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Abstract. In this article we consider a graphene sheet that is folded in various 
compact geometries with arbitrary topology described by a certain genus, g. While 
the Hamiltonian of these systems is defined on a lattice one can take the continuous 
^/-^ \ limit. The obtained Dirac-like Hamiltonian describes well the low energy modes of 

f^**) . the initial system. Starting from first principles we derive an index theorem that 

£SJ ' corresponds to this Hamiltonian. This theorem relates the zero energy modes of 

the graphene sheet with the topology of the compact lattice. For g = and g = 1 
these results coincide with the analytical and numerical studies performed for 
fullerene molecules and carbon nanotubes while for higher values of g they give 
predictions for more complicated molecules. 
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i 1 1 Introduction 

> ■ 

The spectrum of graphene and its various geometrical configurations has been the focus of ex- 
tensive study [1I2I3I4I5I6] . It provides a physical system where a unique interplay is witnessed 
between geometry and electronic properties such as conductivity. Nevertheless, a unified picture 
has not been derived so far due to the richness in behavior of the various geometrical config- 
urations as well as the difficulty in approaching them analytically. One of the interests is to 
study the number of electronic eigenstates with zero energy that determine the conductivity of 
the system and its ground state degeneracy. Previous methods for obtaining the zero modes of 
the system are based on lengthy analytical or numerical procedures. As a possible alternative 
the much celebrated index theorem [7] offers an analytic tool that relates the zero modes of 
elliptic operators with the geometry of the manifold on which these operators are defined. This 
theorem has a dramatic impact on theoretical and applied sciences jSJ. It provides information 
about the spectrum of widely used elliptic operators based on simple geometric considerations 
that could be otherwise hard or even impossible to determine. 

In this article we would like to describe the effect geometrical deformations have on the 
spectrum of graphene. For that we shall establish a version of the index theorem [7|8|9j that 
relates the number of zero modes of graphene wrapped on arbitrary compact surfaces to the 
topology of the surface. In our pursue we shall ignore changes in the couplings caused by 
the geometrical deformations and we shall focus only on the effect the geometry has on the 
spectrum of graphene. As we shall see our results are in good agreement with the known cases 
of icosahedral fullerene molecules [10] and graphite nanotubes [TT] where the spectrum has been 
determined analytically or numerically. Similar approaches for the ground state degeneracy of 
fractional quantum Hall systems in the planar case or on high-genus Riemannian surfaces have 
been taken in [12113114115116] . 
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2 The graphene sheet 

First we shall present an overview of the properties of a flat sheet of graphene. When considering 
its low energy limit a linearization of the energy is possible due to the presence of individual 
Fermi points in the spectrum. This results in a Dirac equation [3] , which describes well the low 
energy behavior of the system. 

Graphene consists of a two dimensional honeycomb lattice where Carbon atoms occupy 
its vertices. When we adopt the tight-binding approximation the model reduces the system of 
coupled fermions on a honeycomb lattice [4] (see Fig. [J) . The relevant Hamiltonian is given by 

H = J E w 

<i,3> 

where J > denotes the tunneling coupling of the electrons along the lattice sites, <i,j> 
denotes nearest neighbors and at, are the fermionic creation and annihilation operators at 
site i with the non-zero anticommutation relation {a^, Oj} — <5y. The original lattice can be 
split into two triangular sublattices, A and B, that correspond to the black and blank circles in 
Fig. [1] This facilitates the evaluation of the dispersion relation of graphene, which is given by 



£(p)=±jyi + 4cos 2 ^+4cos^cos^, (2) 

where the distance between lattice sites is normalized to one. By solving the equation E(p) = 
one deduces that, at half-filling, graphene possesses two independent Fermi points, denoted by 
K_|_ and K_, instead of Fermi lines. This rather unique property makes it possible to linearize 
its energy by expanding it near the conical singularities of the Fermi points. It is not hard to 
show that by restricting near the Fermi points the resulting Hamiltonian takes the form of the 
Dirac operator 

3 T 

H ± = ± T 7> Q , (3) 

where repeated indices are summed over the spatial cooridiantes x, y. The Dirac matrices, 
7°, are given by the Pauli matrices, 7° = a a , and ± corresponds to the two independent 
and oppositely positioned Fermi points. Hence, the low energy limit of graphene is described 
by a free fermion theory. Eigenstates of this Dirac operator are two dimensional vectors, called 
spinors, given by (| K.±A) , | K±B)) T , where A and B denote the two sublattices and K± denote 
two independent Fermi points chosen such that K_ = — K + . 




Fig. 1. The honeycomb lattice comprises of two triangular lattices, A, denoted by black circles and, 
B, denoted by blank circles. A single pentagonal deformation can be introduced by cutting a 7r/3 sector 
and gluing the opposite sites together. 
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3 Curvature deformations and effective gauge fields 

Here we are interested in surfaces with arbitrary topology so we need to introduce curvature in 
the initially flat honeycomb lattice. This is achieved by selectively inserting lattice deformations. 
In doing so, we shall demand that each lattice site has exactly three neighbors and that the 
lattice is inextensional that is it is free to bend, but impossible to stretch. The minimal alteration 
of the honeycomb lattice that can introduce curvature without destroying the cardinality of the 
sites is the insertion of a pentagon or a heptagon; this corresponds to locally inserting positive 
or negative curvature, respectively. Other geometries are also possible, leading to similar results 
as we shall see in the following. 

To introduce a single pentagon in a honeycomb lattice, one can cut a 7r/3 sector and glue 
the opposite sides together, as illustrated in Fig. Q] This causes no other defects in the lattice 
structure. We shall demand that the spinors are smooth along the cut remedied by introducing 
compensating fields which negate the discontinuity |5l6j . Indeed, the cut introduced in Fig. [T] 
causes an exchange between A and B sublattices. This discontinuity can be remedied by intro- 
ducing in the Hamiltonian the non-abelian gauge field A circulation 

j> A^dx" = ^r y (4) 

where r y is the Pauli operator that mixes the K + and K_ spinor components. This flux can 
be attributed to a fictitious magnetic monopole inside the surface with a charge contribution 
of 1/8 for each pentagon [17j . In addition, moving a frame around the pentagonal deformation 
gives a non-trivial coordinate transformation. The effect of this transformation on the spinors 
can be described by a spin connection f2. This is chosen such that its flux around the pentagon 
is given by 

j n^dx* = -~a z (5) 

and measures the angular deficit of 7r/3 around the cone. 

The modified Dirac equation, which incorporates the curvature and the effective gauge field, 
couples the K± spinor components together due to the non-abelian character of A. Since this 
is the only mixing term they can be decoupled by a single rotation that gives 

3 7 

—r{p» - in„ - iA k ^ k = (6) 

where k = 1, 2 denotes the components in the rotated basis with the circulation of the abelian 
now field given by § A^dx^ = ±7r/2 (no summation is considered in k). The curved space 
Dirac matrixes 7 M are given by 7^ = cr Q e^, where is the zweibein of the curved surface 
with metric g^ v that defines the local flat reference frame, 77^,3 = e^e^g M „ . They satisfy the 
anti-commutation relations {7^,7^} = 2g piy , where g^ v is the inverse of g^ v . The curvature of 
the surface is given by the tensor 

where the Christoffel symbols are defined by 

r lv = ^9 ap (d^ P + d u g w - dpg^) 

The Ricci tensor is given by lZp V = lZ a p Va and the scalar curvature is given by 1Z = g^TZ^. 
The field strength that corresponds to the abelian gauge potential, A*, is given by = 
dpA^ — d„A k . Equation ([6]) faithfully describes the low energy behavior of graphene, such as 
its zero modes, when it is deformed to an arbitrary surface. 
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4 Index theorem and graphene 
4.1 The index theorem 

Since the obtained Dirac operator is an elliptic operator, it is possible to employ the index 
theorem 7 8 9| to gain information about its low energy spectrum. Indeed, the index theorem 
gives an insight in the structure of the spectrum of certain operators without the need to 
diagonalize them. This information can be derived from general properties of the operators and 
the geometry of the space, M, they are defined on. A two dimensional Dirac operator defined 
on a surface coupled to a gauge field can be given by the general form 

pt 

P 

where P is an operator that maps from a space V+ to the space V- , while P' maps from V_ to 
V+. As we are interested in the zero modes, we can define the dimension of the null subspace 
of P and P^ by v+ and and i/_ respectively. To facilitate the bookkeeping we introduce the 
chirality operator 75 by 

'l 
-1 



75 



so that it anticommutes with the Dirac operator and it has the states in V± as eigenstates 
with corresponding eigenvalue ±1. As we are interested in the zero modes we can consider the 
operator^) 2 , which is diagonal 

¥ \ PP\ 

and has the same number of zero modes as p. One can easily show that the operators PP^ 
and P^P have the same number of non-zero eigenstates. Indeed, if there is a state u such that 
PP'u = Xu then the state P'u is an eigenstate of the operator P'P with the same eigenvalue, 
PtP(Pt-it) = X(P'u). In order to demonstrate the index theorem we shall employ the heat 
kernel expansion method [18] . Consider a two dimensional compact surface, M. Then one can 
consider the expansion 

Tr(/V* 6 ) = ^-^« fc / 2 a fe (/, J D) (7) 

k>0 

where Tr denotes the trace of matrices and the integration of spatial coordinates. a n (f, D) are 
the expansion coefficients that one needs to determine as a function of the operators / and D. 
For / = 75 and D = p 2 one deduces that 

Tr( 75 e-^ 2 ) = Tr(e~ tptp ) - Tr(e~ tppt ) = J>-* Al ) - ]>>- tA2 ) 

(AO (A a ) 

where Ai and A2 are the eigenvalues of the operators P^P or PPT respectively. But we have 
shown that for every eigenstate of the operator P^P there is a corresponding eigenstate of PP^ 
with exactly the same eigenvalue. Thus only the zero eigenvalues remain giving 

Tr( 75 e-^ 2 ) =i/ + -i/_ 

Combining this result with relation ((TJ it is possible to deduce, that, for / = 75 and D —p 2 , all 
of the coefficients should be zero except for 02 where a 2 = v + — v_ = index(p). The value of 
ci2 defined from ([7]) can be found from the first order term in the t expansion of the exponential. 
Considering that p 2 = — g^ u V u + jil^, r f\^Fia> ~ where V M is the reparametrization 
and gauge covariant derivative, one can easily deduce that 



a 2 = Tr 



75(^,71^-^) 



T 
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where T is the held strength, 1Z is the scalar curvature and the integration runs over the whole 
compact surface M. Thus, we have an analytic way to evaluate the index ofp by 



index(p) = 




The absence of the curvature term in this formula is due to the traceless nature of 75 and it is 
a characteristic of two dimensions. If one can evaluate the integral of the field strength over the 
whole compact surface then the least number of zero modes is determined. It is worth noting 
that for compact surfaces this integral is an integer due to the Dirac quantization condition of 
the monopole charges [P7] . 



4.2 Application to graphene 

Our aim is to evaluate the contribution from the gauge field, in ((SJ) for the particular case of 
a folded sheet of graphene in a compact surface. In a previous section we determined how for 
each lattice deformation a gauge field circulation is introduced. If one could determine the total 
number of deformations for a particular compact geometry of the lattice then we would be able 
to determine the indexilp). At this point we shall assume that curvature is introduced by only 
inserting pentagons and heptagons in the lattice. Interestingly, one can evaluate the number of 
such deformations in a lattice necessary to generate a compact surface by employing the Euler 
characteristic. Indeed, for V, E and F being respectively the number of vertices, edges and 
faces of a lattice defined on a compact surface with genus g, the Euler characteristic, Xi is given 

by 

X = V-E + F = 2(1- g). (9) 

Take the total number of pentagons, hexagons and heptagons in the lattice to be, 715, ne and 777, 
respectively. Then the total number of edges is given by E — (5n 5 +6n 6 + 7ri7)/2 as each polygon 
rii contributes i edges, but each edge is shared by two polygons. Similarly the total number of 
vertices and faces can be evaluated to be V = (5ri5 + 6ne + 7nr)/3 and F = 715 + tiq + n-r, giving 
finally 

n 5 -n 7 = 12(1 -g). (10) 

This result signifies that non-trivial topologies necessarily introduce an imbalance in the num- 
bers of pentagons and heptagons. Moreover, inserting equal numbers of pentagons and hep- 
tagons do not change the topology of the surface as they cancel out. This is consistent with the 
effective gauge field description where a pentagon and a heptagon have opposite flux contribu- 
tions. As particular examples we see that Eqn. (jTOj) reproduces the known case of a sphere with 
.9 = giving x = 2 and a number of defects n 5 = 12 and n 7 = 0. This is the lattice of the Ceo 
fullerene. For the torus we have g = 1 for which x — and 715 = 777 = reproducing the lattice 
of the nanotubes. For a genus-2 surface we have x = — 2, 71,5 = 0, 717 = 12. In all these examples 
equal numbers of pentagons and heptagons can be inserted without changing the topology of 
the surface. 

Now we are in position to evaluate the index(p). The contribution of the gauge field term 
in ([8|) can be calculated straightaway from the Euler characteristic. It is obtained by adding 
up the contributions from the surplus of pentagons or heptagons. Thus, the total flux of the 
effective gauge field can be evaluated by employing Stokes's theorem, giving 

h!l?=h a = 2» s - ^ = ^ - ^ 

where the sign ± corresponds to the k = 1, 2 gauge fields and the summation runs over all the 
surplus of pentagons or heptagons. Hence, from ([5]). one obtains 



6 



Will be inserted by the editor 



Consequently, the least number of zero modes is given by 6|1 — g\, which coincides with their 
exact number if V- = or u + = 0. This is actually what happens in most of the cases when 
the index is non-zero. 



4.3 Zero modes for fullerenes and nanotubes 

The above result relates the number of zero modes of a graphene sheet with the genus of the 
surface it has been folded. As expected it reproduces the number of zero modes for the known 
molecules. The fullerene, for which genus g — 0, has six zero modes which correspond to the 
two triplets of Cgo and of similar larger molecules [4119] . For the case of nanotubes, we consider 
periodic boundary conditions, which give effectively a torus with g = 1. In this case, formula 
(|12| gives v + — V- =0. This is in agreement with previous theoretical and experimental results 
[20111) which have the nanotubes with either no zero modes (e.g. zigzag nanotubes) or with zero 
modes that satisfy v + = V- (e.g. armchair nanotubes). This is a consequence of the symmetry 
between the two opposite directions along the nanotube. 

5 Conclusions 

In this article we have employed the index theorem [7] to enumerate the zero modes of a graphene 
sheet when compactified on arbitrary genus surfaces. Our results are in good agreement with 
the presently studied cases of fullerenes and nanotubes. The only approximation employed here 
was the continuous limit for obtaining the Dirac operator. This approximation is valid if we 
restrict to the low energy spectrum of the system described by large wavelengths. In this limit, 
the lattice spacing or the conical singularities of the pentagonal deformations do not affect the 
low energy modes. Thus, larger fullerene molecules than the Cgo provide more accurate results. 

As an additional example we can consider a graphene sheet folded on an octahedron. In this 
case six square plaquettes have to be inserted in the honeycomb lattice. As square plaquettes do 
not have the frustration properties of pentagons, no effective gauge field is introduced. Thus, 
the index in this case is zero, agreeing with previous considerations fj]. Beyond the known 
examples the version of the index theorem presented here gives a relation between the zero 
modes of more complex molecules. Even if the latter are defined on compact surfaces they are 
related to experimentally relevant ones by imposing appropriate periodic conditions. Thus, the 
metallic properties of these molecules can be induced by simple geometrical considerations. 

This research was supported by the Royal Society. 
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